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ABSTRACT 

I present detailed analysis of the near infrared spectrum of HD 3651B, a faint, 
co-moving wide companion to the nearby planet-hosting star HD 3651. These 
data confirm the companion as a brown dwarf with spectral type T8, consistent 
with the analysis of Luhman et al. Application of the semi-empirical technique of 
Burgasser, Burrows & Kirkpatrick indicates that HD 3651B has T^fj = 790±30 K 
and logs' = 5.0±0.3 for a metallicity of [M/H] = 0.12±0.04, consistent with a 
mass M = 0.033±0.013 Mq and an age of 0.7-4.7 Gyr. The surface gravity, 
mass and age estimates of this source are all highly sensitive to the assumed 
metallicity; however, a supersolar metallicity is deduced by direct comparison 
of spectral models to the observed absolute fluxes. The age of HD 365 IB is 
somewhat better constrained than that of the primary, with estimates for the 
latter ranging over ~2 Gyr to >12 Gyr. As a widely orbiting massive object to 
a known planetary system that could potentially harbor terrestrial planets in its 
habitable zone, HD 365 IB may play the role of Nemesis in this system. 

Subject headings: Galaxy: solar neighborhood — stars: individual (HD 3651B, 
Gliese 570D) — stars: fundamental parameters — stars: low mass, brown dwarfs 
— stars: planetary systems 



1. Introduction 



Luhman et al.l (120061 ) and iMugrauer et al.l (120061 ) have recently reported the discovery 



of a co-moving, widely-separated (43"; 480 AU projected separation) companion to the 
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nearby KO V dwarf HD 3651. Low resolutio n near- infrared spectr oscopy of this source has 
identified it as a late T-type brown dwarf (iLuhman et al.l l2006l ). The HD 3651 system 
also hosts a closely-separated (a = 0.296±0.017 AU), high-eccentricity (e = 0.64±0.04) 
sub-Saturn-mass pl anet (Msini = 0.23 ib0.02 MjupHer), HD 3651b, identified from radial 
velocity variability (IFischer et al.ll2003r). While several exoplanet-host stars are known to 



have s tellar- mass companions (e.g., 



2OO2I : lEggenberger. Udry &: Mayor 



owrance. Kirkpatrick fc Beichmarul2002l : iPatience et al. 



2004 



Raghavan et al.ll2006l ). HD 365lABb is the first 



star/brown dwarf /planet system to be identified, and is therefore a key target for studies of 
all three of these mass-delineated classes. 

In this article, I analyze low-resolution near-infra r ed sp ectroscopy of HD 3651B, us- 
ing the technique of iBurgasser. Burrows fc Kirkpatrickl (120061 . hereafter BBK) to derive the 
effective temperature (Te//), surface gravity (log (7), metalhcity ([M/H]), mass and age of 
this source. Acqui sition and characterization of the spectral data, obtained with SpeX 
faavner et aDl2003h mounted on the 3m NASA Infrared Telescope Facility (hereafter IRTF) 
is described in § 2. Spectral analysis to derive the physical properties of HD 3651B is pre- 
sented in § 3, with the importance of metallicity effects highlighted by the determination 
that this companion must be metal-rich like the system's primary. Finally, in § 5 I speculate 
on the role of HD 3651B as a Nemesis analog in the HD 3651 system. 



Observations 



HD 365 IB was observed with SpeX on 3 September 2006 (UT) under clear and dry 
conditions with good seeing (O'.'T at J-band). The target was acquired using the SpeX guiding 
camera, and low resolution near infrared spectral data were obtained using the SpeX prism 
mode and 0'.'5 slit, yielding a spectral resolution A/AA ~ 120. The slit was aligned to the 
parallactic angl e to mitigate d i fferent ial refraction effects. Note that this setup was similar 
to that used by iLuhman et al.l (120061 ) for their spectral observations of HD 3651B, with the 
exception of a slightly narrower slit (0'.'5 versus O'.'S). Eight exposures of 180 s each were 
obtained in an ABBA dither pattern along the slit. The AO V star HD 7215 was observed 
immediately afterward and at a similar airmass (1.06) for flux calibration. Internal flat fleld 
and Ar arc lamps exposures were also acquired for pi xel response and wavelength calib ration. 
All data were reduce d using SpeXtool vers ion 3.3 (ICushing. Vacca. fc Rayned 120041 ) using 
standard settings (cf. IBurgasser et al.ll2004l ). 



Figure [T] compares the reduced sp ectrum of HD 365 IB to equivalent data for the T7.5 
companion brown dwarf Gli ese 570D (jBur gasser et al.ll2000l ) and the T8 fleld brown dwarf 
2MASS J04151954-0935066 JBurgasser et al..i200a . hereafter 2MASS 0415-0935). The strong 
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H2O and CH4 bands pres e nt in the spectrum of HD 3651B are consistent witli tlie obser- 
vations of iLufiman et al.l (120061 ) . and confirm tliis object as a la te T-type brown dwar f. 
Classification of this source was done using the spectral indices of IBurgasser et al.l (120061 ). 
The measurements and associated subtypes are as follows: H2O-J = 0.044 (T8), CH4-J 
= 0.203 (T7.5), H2O-H = 0.186 (T8), CH4-H = 0.134 (T7.5) and CH4-K = 0.043 (T8), 
for a mean cla s sificat ion of T8. This value is consistent with the T7.5 subtype derived by 
Luhman et al. J2006h given the typ ical 0.5 subtype uncertainty in near infrared T dwarf 



classifications (IBurgasser et al.ll2006r ). 



3. The Physical Properties of HD 3651B 
3.1. Derived Properties and Metallicity Dependence 

As a companion to a well-studied, nearby star, the properties of H D 365 IB could be 



readi l y derived by adopt i ng the age and meta llicity of the primary (e.g. iKirkpatrick et al. 



2OOII : IWilson et al.l I2OOII : iLuhman et al.l |2006| ) . Alternately, the semi-empirical technique 
described in BBK can be employed. This method involves the comparison of H2O and 
color spectral ratios measured on the spectrum of a late-type T dwarf, which are separately 
sensitive to Teff and logg (for a given metallicity), to the same ratios measured on theoretical 
spectral models. The latter are calibrated to r eproduce the measured indices for the near 
infrared SpeX prism spectrum of Gliese 570D (IBurgasser et al.ll2004l ). which has reported 
physical parameters of T^// = 782-821 K, log c/ = 4. 95-5.23 and [Fe/H] = 0.09±0.04, based o n 
empirical measurements and evolutionary models (iGeballe et al.l 1200 ll : ISaumon et al.ll2006l ). 
This method provides constraints on both Teff and logg which can then be used to infer 
mass and age with evolutionary models or empirically via the source's bolometric luminosity. 



In BBK and iLiebert fc Burgasseii (120061 ). this method was applied assuming solar or 
subsolar metallicities. Ho wever, the primary of the HP 36 51 has a metallicity greater than 
sola r, [Fe/H] = 0.12±0.04 dSantos. Israehan fc Mayorll2004j ). It has been shown in both BBK 



and lLiebert &: Burgasserl (120061 ) that slightly subsolar metallicities can significantly skew the 
derived Teff and logg values. Hence, several cases were con sidered for < [M/H] < 0.24. 
The spectral models of iBurrows. Sudarsky &: Hubenyl (120061 ) were employed, spanning 700 
< Te// < 1100 K, 4.5 < logg < 5.5 (cgs) and < [M/H] < 0.5. Four pairs of indices were 
compared between (H2O-J, H2O-H) and (K/J, K/H)Q 

Figure [2] illustrates the {Teff, logg) phase spaces for both [M/H] = and 0.12 over which 



^For definitions of these spectral ratios, see BBK and IBurgasser et al.l (|2006l ). 
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the spectral ratios H2O-J and K/H measured on the spectrum of HD 3651B match cahbrated 
values from the models. Their intersection yields robust constraints on these parameters for 
the source. All four pairs of indices result in a combined Teff = 790±30 K and logg = 
5.0±0.2 for [M/H] = 0.12 (Table 1); including the full uncertainty in [M/H] increases the 
uncertainty on logg to 0.3 dex|^ For [M/H] = 0, the analysis yields Teff = 810±30 K and 
logg = 4.7±0.2. Note that a decrease in the assumed metallicity leads to a decrease in 
the der ived surface gravity. This is consistent with trends for subsolar metallicity T dwarfs 
(BBK; iLiebert fc Burgassen 120061 ) . and is largely due to the impact of both metallicity and 
surface gravity on the relative opacity of c ollision-induced H2 absorption that dominates the 



i^'-band opacity (e.g.. ISaumon et al.lll994l ). 







The T^ff and log^ values of HD 3651B for [M/H] = 0.12±0.04 infer M = 0.03 3 ±0.01 3 M, 
and an age r = 0.7-4.7 Gyr, based on the evolutionary models of iBurrows et al.l (I2OOII ). For 
solar metallicity, these values decrease to M = 0.020±0.005 M© and r = 0.4-1.2 Gyr, primar- 
ily due to the lower surface gravity indicated in this case. An estimate of mass can also be 
obtained from the bolometric lu minosity (L) of HD 36 51B by the relation M = Lg/AirGcrT^jj-. 



Using the luminosity derived by lLuhman et al.l ((20061), logL/L© = — 5.60±0.05, and the Teff 



and logg values derived from the spectroscopic analysis yields M = 0.023^q;q^q Mq (M = 
0.01ll^:J^|]^ M0) for [M/H] = 0.12 ([M/H] = 0.0). These values are consistent with those 
derived from the evolutionary models, albeit with larger uncertainties driven primarily by 
the uncertainty in logg. 



3.2. HD 3651B: A Metal-Rich T Dwarf 

The method outlined above, based on only two indices, cannot discriminate between 
[M/H] = and 0.12 for HD 3651B, for in both cases it is possible to derive robust con- 
straints on both Te// and logg (this is not always the case; cf. Figure 8 of BBK). A third 
spectral index that is uniquely sensitive to metallicity is required, and BBK have suggested 
the ratio between the 1.05 and 1.25 /im peak fluxes [Y/J index) as a possible option. 
This ratio is sensitive to the far red wings of the pressure- broadened 0.77 /im K I dou- 
blet, and models indicate considerab le metallicity sensitivity to this feature. However, the 



Burrows. Sudarsky fc Hubenyl (120061 ) models do not reproduce the observed trends between 



F-band and J-band, likely due to current uncertainties in the treatment of the far-wing 



^Note that additional systematic uncertainties of ±50 K and ±0.1 dex in the derived Tg// and logg 
values may be present due to uncertainty in the parameters of Gliese 570D and possible systematic effects 
in the spectral models (BBK). 
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opacities (IBurrows fc Volobuyevll2003l : lAUard et al.ll2003l ) 



An alternate procedure, used by lLiebert &: Burgasserl (120061 ) . is to compare the observed 
absolute fluxes of the brown dwarf to spectral models based on the best fit parameters for a 
given metallicity. This is done in Figure [3l which s hows the absolute flux -calibrated spectrum 
of HD 3651B, based on J-band photometry fromlLuhman et"aL ( 2006 . 16.16±0.03) and the 
parallax of HD 3651A from HIPPARCOS JPerrvman et all IiQQtI . tt = 90.0±0.7 mas); to 
predicted fluxes for two Tg// = 800 K models, one with logg = 4.7 and [M/H] = and one 
with logg = 5.0 and [M/H] = 0.12. While the models provide poor fits to the 1.1 and 1.6 fim 
CH4 bands, the result of well-documented deficiencies in CH4 opacities at these wavelengths, 
it is clear that the lower surface gravity, solar metallicity model is too bright as compared to 
the data. The [M/H] = 0.12 model, on the other hand, provides a reasonably good match 
to the absolute fluxes at each of the spectral peaks. This comparison clearly supports a 
supersolar metallicity for HD 3651B. 

Another consideration is the inferred age of HD 3651B as compared to previ ous age de- 
termin ations for HD 3651A. Published estimates for the primary vary appreciably. IWright et al 
(I2OO4J ) report an age of 5.9 Gyr on the bas i s of w eak Ca II HK emission (logR'fff^ = —5.07) 
and t he age/activity relation of iDonahud (119931 ). On the other hand 



Rocha-Pinto et al. 



from 
has 



(12004 ) report a chromo s pheri c age of only 2.1 Gyr based on a higher value of log R'hk — ~4.85 



Soderblom et al.l (119931 ). Indeed, a 35% variation in Ca II HK emission of HD 3651B 



3een seen over 17 years of observations at Mount Wilson Observatory (IDuncan et al. 



19911). From isochrone fi tting, iNordstrom et al.l (l2004f) derive a minim um age of 2.6 Gyr, 
Valenti &: Fischerl ( 20051 ) an age of 8^^ Gyr and iTakeda et al.l (12006h a minimum age oi' 
11.8 Gyr. While there appears to be little consensus, an older age for HD 36 51A is most 



consi stent with its chromospheric i nactivity, its long rotation period (48 days; iNoyes et al. 



19841 ) and its UVW space velocities (lEggenlll989l ). A value in the range 2-12 Gyr is indicated. 



This result again argues in favor of HD 365 IB being metal-rich. The derived age for 
the solar metallicity case is considerably younger than that of the primary; while for [M/H] 
= 0.12±0.04, the ages of the primary and secondary overlap in the range of 2-5 Gyr. Note 
that this age constraint is more precise than that inferred from the primary alone, and is 
consistent with both components having been formed at the same time and from the same 
gas/dust reservior (same age and metallicity) 

There are two conclusions that may be drawn from this analysis. First, it is clear 
that the determination of Tg// and logg, and hence mass and age, of an individual brown 
dwarf using the BBK technique is highly sensitive to the assumed metallicity, even for small 
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variationsjfl This reflects the overall sensitivity of cool brown dwarf spectra to changes in pho- 
tospheric abundances, and is not surprising given that the entire near infrared spectrum of a 
T dwarf is draped by atomic and molecular absorptions. While metallicity effects have been 
previously noted in the spectra of metal-p oor T dwarfs such as 2MASS J0937348 7+2931409 



(Burgasser et al. 



1999 



2002 



Knapp et al.ll2004 BBK) and 2MASS J12373919+6526148 (IBurgasser et al. 



Liebert fc Burgasserl 120061 ). it is now clear that metal-rich T dwarf spectra can also 



exhibit metallicity effects. These results emphasize the importance of characterizing the 
overall abundances — and abundance patterns — in deriving a complete and accurate char- 
acterization of a low-temperature brown dwarf. Note that HD 3651B and Gliese 570D may 
be a particularly useful pair in studying metallicity effects, as they appear to be identical in 
nearly all aspects {Teff, log (7, age, mass and spectral energy distribution) except metallicity. 

Second, it is possible that coeval brown dwarf companions may provide additional or 
even better constraints on the ages of their host systems than some stellar primaries. While 
age estimates for stars on the basis of chromospheric activity and kinematics can vary ap- 
preciably, and are less discriminating at later ages, the example of HD 365 IB demonstrates 
that a late-type brown dwarf can potentially be age-dated with higher precision. However, 
the accuracy of these ages remains unclear, as they are tied to both the fidelity of evolution- 
ary models and, in the case of the BBK technique, the accuracy and calibration of spectral 
models. In the former case, there have been few empirical checks on the evolutionary models 
of brown dwarfs at late ages0 although BBK found consistency between the radii of field T 
dwarfs empirically derived from luminosity measurements and those based on evolutionary 
models. In the latter case, there appears to be consistency between different sets of spectral 
models, as an equi valent analysis of the spectrum of HD 3651B using the COND models of 



AUard et al.l (1200 ll ) produced nearly identical results. Hence, at first glance it appears that 



late-type brown dwarf companions could be useful sources for characterizing the properties 
of a presumably coeval system, although further empirical tests of brown dwarf evolutionary 
models are needed. 



■^[M/H] = 0.12 implies only a 32% increase in metal abundance over solar. 

^At very young ages (~l-5 Myr), both iMohantv et al.l |2004h and Istassun. Mathieu k Valentil |2006[ ) 
have found discrepancies between evolutionary models and observed parameters. It is possible that these 
discrepancies ar ise from the initial c onditions assumed in the evolutionary models, and may be unimportant 
at late ages (cf. iMarlev et all 120061 ). 
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HD 3651B: An Extrasolar Nemesis? 



HD 365 IB is tlie first wide brown dwarf companion found in a system hosting at least 
one known planet. Tfiis raises an intriguing question: liow lias the evolution of the HD 3651 
planetary system — and specifically the development of life on any putative habitable plan- 
ets in this system — been influenced by this distant brown dwarf companion? This spec- 
ulative question is an appro priate one, as putative terrestrial bo dies within a portion of 
HD 365lA's habitable zone (IKasting. Whitmire fc Revnoldd 119931) can persist without dy- 
namical disruption by HD 3651b ( Jones. Underwood fc Sleep|"2005), e ven considering this 
giant planet's probable inward migration ( Mandell fc Sigurdsson " l2003 ). Terrestrial planet 
formation around the primary of a wide b inary system has also been shown to be feasible 
(IHolman. fc Wiegertl 119971 : iQuintanal |2004| ) . However, the formation, existence and habit- 
ability of terrestrial planets in the HD 3651 system may be inhibite d by dynamical pertur - 
bations from HD 3651B, both directly through orbital modulatior@ fiTakeda fc Rasiol 120051 ): 
or the scattering of small bodies, resulting in sustained impacting onto the planet's surface. 
Yet terrestrial planets on mildly ecc entric orbits can niaintai n habitability depending on 
the properties of their atmospheres (lAdams fc Laughlinl |2006| ) . while increased impacting 
rates may in fact facilitate the emergence of life. In the case of the Earth, the late heavy 
bombardment period between 4.5 and 3.8 Gyr ago resu lted in one particularly cataclysmic 
impact that formed the Moon (IHartmann fc Davisll 19751 ). whose tidal forces have maintained 
Earth's tilt to the Sun and have thereby reduced climate variation. Inipacts of icy bod ies 
may have brought necessary water ice to Earth's surface ( Chybal 119871 : iDelsemmd Il998l ) as 
well as chemical precursors to biotic life Jchvba et al.lll99oh . On the o ther hand, later im - 



pacts clearly have a ne gative effect on the evolution of macrobiotic life (jAlvarez et al 
Hildebrand et al.lll99lh . 



1980 



As an instigator of planetary impacts, HD 365 IB may pla y the role of Nemesis in the 



HD 3651 system. 



his h ypothesized companion to the Sun (IDavis. Hut fc Mullerl Il984 



Whitmire fc Jacksonlll9841) has been proposed to exp 



extinctions (IFischer fc Arthur 



1977 



tering rates (I Alvarez fc MuUei 



1984 



ain apparent periodi cities in massive 



Raup fc Sepkoskilll984: iHut et al.lll987r) . terrestrial cra- 
Rampino &: StothersI 1984 ) , reversals of Earth's mag- 



netic fields (|Raudll985[) and, more recently, the peculiar orbits of inner Oort Cloud plane- 
toids (IBrown. Trujillo fc Rabinowitzll2004l : iGaudi fc Bloonj|2005[ ) . If HD 3651B were present 
at its observed separation but bound to the Sun, it would be invisible t o the naked eye 
{My ~ 27 implies ~ 9) but readily detectable by Hipparcos, Tycho (jPerryman et al. 



^Indeed, Luhman et al. 20061 suggest that the eccentricity of HD 3651b may be driven by secular pertur- 
bations from HD 365 IB. 
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19971) and early wide-field near-inf rared sky surveys {K ~ —1; cf. the Two-Micron Sky Sur- 
vey, Neugebauer fc Leighton 79691) . If HP 3651B was a t the distance from the Sun proposed 
for Nemesis (~1.5xl0^ AU; iDavis. Hut &: Mullerlll984j ). it would be invisible on visual pho- 
tographic plate images (V ~ 20) but a relatively bright source (J ~ 10) in the Two Micron 



All Sky Survey fISkrutskie et al.ll2006l ). 



Searches for a solar Nemesis have so far turned up negative (e.g., iPerlmutterl Il986l ). 
However, as widely-separate d low mass stellar corn panions appear to be prevalent in known 
planet-hosting systems (e.g.. iRaghavan et al.ll2006l ). and with widely-separated brown dwarf 
companions to these systems only now being identified. Nemesis counterparts may be a 
common facet of all planetary systems and an important consideration for the emergence of 
terrestrial-based life. 



The author thanks D. Chakrabarty and J. Winn for their helpful comments on an early 
version of the manuscript, and acknowledges the support of telescope operator E. Volquard- 
sen and instrument specialist J. Rayner during the IRTF observations. The author also 
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Table 1. Derived Physical Properties of HD 3651B. 



Parameter 


[M/H] = 0.00 


[M/H] = 0.08 


[M/H] = 0.12^ 


[M/H] = 0.16 


[M/H] = 0.24 


Te// (K) 


810±30 


790±30 


790±30 


790±30 


780±40 


log 5 (cgs) 


4.7±0.2 


4.9±0.2 


5.0±0.2 


5.0±0.2 


5.2±0.2 


Mass (Mq) 


0.020±0.005 


0.027±0.007 


0.032±0.008 


0.035±0.011 


0.047±0.013 


Age (Gyr) 


0.4-1.2 


0.7-2.8 


1.0-3.7 


1.0-4.7 


2.6-9.3 



Santos. Israelian &: Mayor 



(|200J) derive a metallicity for the primary of [Fe/H] = 0.12±0.04. 
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T 




Wavelength (|im) 

Fig. 1. — Prom top to bottom, normalized near- infrared spectra of the T7.5 companion 
Glicse 570D, HD 3651B and the T8 field brown dwarf 2MASS 0415-0935, all obtained with 
SpeX on IRTF. All spectra are normalized at 1.27 fim and offset by constants (dotted 
lines). Major spectral features are noted; the © symbols designate regions of strong telluric 
absorption. 
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Fig. 2. — Phase spaces of Tg// and \ogg constrained for HD 3651B by the indices H2O-J 
and K/H using the technique of BBK for metaUi cities [M/H] = 0.0 (left) and . 12 (ri ght). 
Model indices are measured from calculations by iBurrows. Sudarsky &: Hubenyl (120061 ) and 
calibrated to observations of Gliese 570D. The intersecting cross-hatched regions provide 
robust constraints on Teff and logg. Note the higher gravity constraint derived for the 
higher metallicity (see also Table 1). 
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Fig. 3. — Comparison of the absolute flux- calibrated spectrum of HP 36 5 IB (b lack line) 
to two Te// = 800 K spectral models from iBurrows. Sudarsky fc Hubenyl (120061 ): logg = 
4.7 and [M/H] = 0.0 (red dashed line), and log^ = 5.2 and [M/H] = 0.12 (blue triple-dot 
dash ed line). Flux calibrat ion of the HD 365 IB spectrum is based on J-band photometry 
from iLuhman et al.l (120061. 16.16=fc0 .03) and the parallax of the primary as measured by 
HIPPARCOS (jPerryman et al.l 119971 . tt = 90.0±0.7 mas). Uncertainty in this calibration at 
the J-band peak is indicated by the error bar. Note the agreement between the data and 
the [M/H] = 0.12 model. 



